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In this work we investigate the inverse of the celebrated Bohigas-Giannoni-Schmit conjecture. 
Using two inversion methods we compute a one-dimensional potential whose lowest N eigenvalues 
obey random matrix statistics. Our numerical results indicate that in the asymptotic limit, TV —¥ 
oo, the solution is nowhere differentiable and most probably nowhere continuous. Thus such a 
counterexample does not exist. 
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I. INTRODUCTION 

One of the main milestones of quantum chaos is the 
celebrated Bohigas-Giannoni-Schmit (BGS) conjecture 
[l|. This conjecture forms a link between the statisti- 
cal properties of the spectra of quantum systems with 
chaotic classical limit and random matrix theory (RMT): 
"Spectra of time-reversal-invariant systems whose clas- 
sical analogs are K-systems show the same fluctuation 
properties as predicted by GOE Jj|/." 

This conjecture establishes a connection between 
chaotic systems and RMT, i.e. Chaos — > RMT. A num- 
ber of numerical evidences corroborate this conjecture 
3, but a complete analytical proof is still missing. As 
pointed out in Ref. Q this universality is a consequence 
of the properties of the long periodic orbits whose semi- 
classical limit is also universal. The intimate relation 
between RMT statistics and its breakdown over energy 
scales of O(hfT) where h is Planck's constant and T is 
the period of the shortest periodic orbit, was emphasized 
in Refs. [1, H[. The first attempt to prove the BGS con- 
jectures is due to Andreev et al. @, Q but later it turned 
out to be incomplete. More recently the link between 
chaotic systems and RMT has been investigated using an 
improved semiclassical technique fc^-fl(| . Finally we have 
to mention that certain chaotic models show Poisson- 
statistics [TTJ] allowing for the possibility of exceptions, 
as well. 

Hence it is indeed very difficult to understand the ex- 
act relationship between quantum chaos and RMT. Cer- 
tainly it should be helpful to investigate the inverse BGS 
conjecture, i.e. does the presence of GOE level statis- 
tics ever lead to chaos in the classical limit? This would 

? 

be the opposite link: RMT— ^Chaos. This question has 
been raised first by Wu et al. [l^ . The authors' aim 
was to find a counterexample, i.e. a one-dimensional 
quantum system with GOE level statistics. However, 
a one-dimensional system is always integrable guaran- 
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teed by Liouville's theorem [Tj]. The choice of d = 1 
is merely because it is easier to do computation in one 
spatial dimension as compared to d = 2 or d = 3. In 
addition the average level separation is of the order of 
0{h d ) which allows a for semiclassical regime in d > 1 
only. Despite the simplicity of d = 1 the task is not at 
all easy, finding a potential with a given spectrum leads 
to an inverse Schrodinger problem. There is almost no 
chance to achieve this analytically, therefore Wu et al. 
used a numerical method to find a one-dimensional po- 
tential with the first ./V energy levels being predefined. 
Their process gave accurate results for the case of GOE 
level statistics, hence they concluded that there is an in- 
tegrable system with GOE level statistics, and provided a 
counterexample to the inverse BGS conjecture. Wu and 
Sprung later carried out a similar calculation [Tij for a 
particular, asymptotically CUE spectrum represented by 
the non-trivial zeros of the Riemann zeta function (see 
also fl5l4l9l |). Their findings led to the concept of a frac- 
tal potential, which may produce a chaotic spectrum. 

Contrary to Wu et al. we can only state that the first 
./V eigenvalues of the potential obtained numerically show 
random matrix type fluctuations leaving any other type 
of behavior for the remaining infinite number of eigen- 
values. Therefore we re- investigate this problem more 
carefully with the emphasis on the N — > oo, limit. 

We examine certain statistical properties of the result- 
ing potential and conclude that no classical system can be 
associated with the constructed quantum potential. Our 
analysis not only confirms a kind of fractal nature of the 
potential, already conjectured by Wu and Sprung (l2l . [l6j . 
but exceeds it by showing that the potential asymptoti- 
cally loses its continuity and differentiability and resem- 
bles the sum of randomly positioned Dirac-dclta func- 
tions. 

The structure of our article reads as follows. In Section 
[II] we describe the numerical methods applied, Section Hill 
compares the results provided by the numerical methods. 
In Section ITVl we analyse the N — >• oo limit and Section 
fVlis left for conclusions. 
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II. THE NUMERICAL METHODS AND 
THEIR PROPERTIES 

We discuss the inverse BGS conjecture and start by cal- 
culating all eigenvalues of a sufficiently large, real, sym- 
metric matrix, thereby creating a GOE spectrum. This 
raw spectrum has to be unfolded in order to guarantee 
uniformity [20| . As a result of this procedure, the average 
distance between two neighboring levels in the unfolded 
spectrum is unity. The first N members of this unfolded 
spectrum are thereafter denoted by e l5 • ■ • , e^r and re- 
main unchanged. 

We apply two techniques, the numerical method used 
als o by Wu et al. , and the dressing-transformation [TBI . 
[Hi H2 m order to find a potential, whose first N levels 
are as close as possible to ei, • • • , e^. By changing N we 
tried to conclude about the N — > oo limit, since in that 
case the entire spectrum should show GOE properties. 
Naturally, infinitely many potentials may have the same 
first N levels, nevertheless their shapes coincide for the 
energy range determined by the prescribed spectrum and 
differ only at higher energies. 

Below we describe the numerical techniques; an op- 
timisation and an iterative one. These two approaches 
provide qualitatively and quantitatively similar results 
concerning the most important properties of the poten- 
tial. 

At first we implemented the same variational algorithm 
introduced by Wu et al. 12j. This method starts with 



an initial guess Vq(x) for the yet unknown potential V(x) 
and one calculates the energy eigenvalues of Vo(x), which 
are denoted by e±, £n- This spectrum may differ 
from the pre-defined energy levels, therefore in the next 
step the potential is altered in a way, that its energy 
eigenvalues fit the expected spectrum, ei, • ■ ■ , ejv, better. 
The same procedure can be continued until convergence 
is reached. 

This method requires an objective function, G, mea- 
suring the distance between the actual and the pre- 
defined spectrum. One may choose G to be 



A' 



G[V(x)]=J2(ei 



(1) 



G implicitly depends on V(x) through the set of eigenval- 
ues, £,;, therefore the task is to minimize G through the 
variation of V(x). Moreover, G is positive semidefinite, 
and equals zero if and only if Ei = ei (i =1,2, . . . , N). 

Employing a standard gradient method the potential 
has to be updated as 



iV 
i=l 



Wz) - V old (x) - 2r, V ( £i - ei ) [^(x)}. (2) 



If the initial potential, Vq(x), is well chosen, the algo- 
rithm converges and provides us with a potential, whose 
first N eigenvalues are very close to the pre-defined ones, 



The average distance in the unfolded spectrum of a 
large GOE matrix is unity, which resembles a harmonic 
oscillator, therefore a parabolic form seems to be a good 
choice for Vq(x). 

Our other method is based on an interesting property 
of the Schrodinger equation: a bare free quantum sys- 
tem can be "dressed up" in an iterative manner, with a 
series of potentials incorporating more and more bound 
states into the system. Thereby this method is generally 
called dressing transformation (l5l [l8l . l2~l| . Eventually 
this algorithm provides a potential, whose first N lowest 
energy eigenstates are bound states with the pre-defined 
energy levels, while all other states are scattering states, 
i.e. ei < e 2 < ■ • ■ < ejy < 0, and = for k > N. 
One may interpret this method differently, and say that 
it shifts N scattering states to negative energies, convert- 
ing them into bound states, while preserving the original 
difference between the levels. It is important to mention 
here, that theoretically this method is exact, and only 
the implementation may introduce numerical error into 
the final solution. 

In the remaining part of this section we comment on 
the uniqueness of the potential, on the speed of the al- 
gorithms, and finally on the numerical error accepted in 
our calculations. 

Uniqueness: If one intends to explore the interaction 
in a quantum problem by analyzing scattering data, it 
is a crucial question whether the experimental scatter- 
ing data is compatible with many different potentials, 
or just with one. For a generic non-relativistic, one- 
dimensional quantum system Borg has rigorously proven 
that one spectrum (eigenvalues for the bound states and 
reflection coefficient for the scattering states) is not suf- 
ficient to uniquely reconstruct the corresponding poten- 
tial (22)]. However, from two spectra (obtained using two 
different boundary conditions) one can determine the po- 
tential uniquely. Restricting ourselves to even potentials 
only, V(x) = V(—x), one spectrum becomes sufficient 
to uniquely determine the potential (23|. Unfortunately, 
noone can use an infinite set of spectrum in a numer- 
ical algorithm, therefore uniqueness is lost apparently. 
Our physical intuition suggests though, that if two po- 
tentials support the same bound states, then these poten- 
tials ought to be identical at least around their bottoms. 

This expectation is valid and justified numerically in 
our calculation. The potential generated by the Wu- 
method and by the dressing transform - after an in- 
significant energy shift - differ only at their "semiclas- 
sical edges" which can be characterized by i t um, where 
V(xtum) = cat, i.e. where a classical particle with energy 
ejy would turn back in the given potential. 

The difference between the potentials is expected out- 
side the [— Xtum, ^turn] interval, since the Wu-method 
should reproduce the parabolic potential characteristic of 
the harmonic oscillator, while in this form of the dressing 
transformation the first derivative of the potential must 
approach zero asymptotically, therefore the potential be- 
comes constant. 
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Speed: According to our experience the Wu method 
was very accurate; the error decreases rapidly in each 
consecutive iteration. Despite its precision, it becomes 
slow for increasing TV mainly due to the repeated calcu- 
lation of every eigenstates, $j(x) (i = 1, 2, . . . , TV). This 
method seems to be prohibitively impractical for TV being 
higher than few hundreds. On the contrary the dressing 
transformation is very fast, but for large TV, its accuracy 
decreased considerably. 

Numerical error: Both methods introduce a numerical 
error due to the discretization. It is, therefore, necessary 
to check whether the potential calculated reproduces the 
spectrum prescribed. If the spectrum e% does not match 
the original set of eigenvalues, e», the appropriate param- 
eters of the algorithms must be adjusted: the variational 
step size, T], in case of the Wu-method, and the discretiza- 
tion step-size of the Rungc-Kutta method for the dressing 
transformation. 

We stopped the algorithms when the average error on 
each energy level has been reduced below a fixed thresh- 
old, namely 

1 N 

-5> t - £l |<l(r 5 . (3) 

i=l 

This seemingly acceptable order of numerical error is eas- 
ily met by the Wu-method, but represents a stringent 
condition for the dressing transform, the later of which 
needed small step-size, ~ 2 x 10~ 4 , for relatively high 
values of TV. We have also investigated the accumulation 
of error and have found that the sum in Eq. © collects 
the main contributions at high energies, i.e. where the 
two methods are expected to deviate from each other. 



III. RESULTS 

This section provides the results. First let us take a 
look at Fig. Q] which depicts the potentials obtained by 
the Wu-method and the dressing transformation for TV = 
50. 

Note that the two methods do indeed give very similar, 
symmetric potentials, as we expected based on Borg's 
theorem. The dressing transformation provides a sym- 
metric potential a priori, while the Wu method is started 
from a symmetric potential which property is inherited in 
each iteration. More importantly, one may separate the 
potential into two parts; a parabola and an oscillatory 
component, the latter of which will be called fluctuation. 

V(x) = Vjx) + Vfi(x). (4) 

Due to the GOE level fluctuations of the first TV eigen- 
values the parabola becomes oscillatory in a finite region 
that is connected to the TVth level, i.e. to the turning 
point Xtum for which V(a;turn) = E. In quantum me- 
chanics the wave function approaches zero very rapidly 
for | a; | > Xtum, thus the oscillatory region is simply the 
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FIG. 1. The potentials calculated with the Wu-method (solid 
line) and the dressing transformation (dashed line) for TV = 50 
are shown. The potentials provided by these methods are 
almost identical in the [-10, 10] range. Inset shows the same 
potentials in the [0, 3] interval. 



interval of [— Xturni ^turo] if E — £jv- Hence the aver- 
age potential becomes wavy in the region according to 
Xtmn corresponding to the TV-th level. By increasing TV, 
the turning point, itum, increases too, therefore in the 
TV — > oo limit the whole potential is expected to be wildly 
oscillatory. This expectation is supported by Fig. [2] al- 
though the largest value of TV used is only 250. However, 
not only the width of the oscillatory region increases, but 
the amplitude of the fluctuation term seems to be larger 
and larger, see Figs. [2]and[3][a). In the following we wish 
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FIG. 2. The potential calculated with the Wu method for 
TV = 5 (bottom), 50 (middle) and 250 (top). The graphs for 
TV = 50 and TV = 250 are shifted vertically. 

to classify the oscillatory component, Vn(x), by measur- 
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ing the density of the fluctuations. For this reason the 
distance between two consecutive extrema could serve as 
a good measure, see Fig. E^b), therefore, we defined Ax 
and Ay as the distances between two adjacent extrema 
in the x and y directions. For a fixed value of TV, the av- 
erage value of Ax and Ay characterizes the fluctuation of 
V(x). The average of Ax is plotted in Figure @] for both 
methods. The graphs suggest a power-law dependence 
on TV. The least-squares fit, using Ax = A + BN~ C 
is also plotted to guide the eye. Provided this depen- 
dence is suitable also in the asymptotic limit, TV — > oo, 
the average width of oscillation in V-n(x) approaches zero, 
since i«0, while the average amplitude of fluctuation, 
Ay, remains positive, see Fig. This asymptotic be- 
havior naturally poses the question: does the potential 
remain differcntiablc as TV — > oo? In order to answer this 
question we introduce a new quantity with the follow- 
ing definition: c = Ay /Ax. This fraction, at any finite 
TV, is not equal to the usual difference fraction, because 
neither Ace, nor Ay is allowed to change continuously as 
in calculus, but here they have a finite and fixed local 
value. For any potential calculated for a fixed TV, we 
have a finite set of c values to analyze statistically. In 
Figure [6] we show the histogram of c for different values 
of TV. One can see that this histogram broadens as TV 
increases. Conclusively, c, the average of c, diverges pro- 
portional to TV, which is seen in Fig. [7J The behavior 
of Ax and c described above allows us to conclude: in 
the TV — > oo limit the whole potential will be rapidly os- 
cillatory, because x turn is growing with TV. Considering 
any finite interval of x, Ax — » as TV — » oo, therefore 
c becomes the derivative of V(x). From the numerical 
observations mentioned above one can conclude (i) there 
are infinitely many extrema, thus infinitely many values 
of c in a finite interval of x, and (m) a finite part of the 
c-s, still infinitely many, diverge with increasing TV, since 
c diverges as well. These findings suggest that the po- 
tential is nowhere differcntiablc. Such a potential cannot 
be compatible with the Hamiltonian formalism of clas- 
sical mechanics, especially with the Hamilton-equation 
p = —dH/dq. Therefore Wu's counterexample for the 
inverse BGS conjecture is not acceptable, since the classi- 
cal limit of a potential whose entire spectrum obeys GOE 
level statistics makes no physical sense. Subsequently, no 
meaning can be assigned to the question whether it leads 
to chaotic dynamics or not. 



IV. MATHEMATICAL PROPERTIES OF THE 
TV -> oo LIMIT 



In order to answer these questions we are going to fo- 
cus on the fluctuation part of the potential and calculate 
for instance its numerical Fourier transform, .F[Vfl(x)]. 
The Fourier components, for different values of TV, are de- 
picted in Fig. [8] It is apparent that the Fourier-transform 
differs from zero only in an interval [0, / max ]- This upper 
bound can be thought of as a characteristic frequency, 
which translates into a characteristic width. As TV in- 
creases / max increases as well (see inset of Fig. [8]) there- 
fore waves with smaller and smaller wavelength appear 
in the graph of the potential. In the TV — » oo limit we 
expect that every wavelength will appear, so the graph 
of the potential is going to lose its characteristic length 
scale. This is exactly a particular property of fractals. In 
Fig. [5] the fractal dimension of the potential is plotted, 
obtained by using the standard box-counting algorithm. 
Both methods employed yield similar behavior: a slowly 
increasing box-dimension as TV increases. Unfortunately 
the computed dimension fluctuates too much to derive 
a definite conclusion. However, the dimension is most 
likely not smaller than 1.6 in the TV — > oo limit. This is 
comparable to the value of 1.5 obtained in an earlier work 
for the potential derived from the nontrivial zeros of the 
Riemann zeta- function [l4j . This means that the poten- 



tial strongly shades the complete plane for large-enough 
TV, which can be seen in the inset of Fig. [5] Earlier we 
concluded that our potential is nowhere differentiable. 
Furthermore continuity would require Ay to decrease and 
approach zero as the TV — > oo. In Fig. [5] one can see that 
Ay does not tend to zero. If Ay increases indefinitely, 
the graph of the potential is not even bounded. In or- 
der to analyse this possibility, we also calculated the area 
under the graph of the potential defined as the definite 
integral of its absolute value / |Vfl(x)| dx. Fig. ITOl shows 
the number of nodes, U, while the inset depicts the total 
area under the graph of the potential of the same graph. 
It is apparent from Fig. [TO] that both the area and the 
number of nodes increases linearly with TV. The area un- 
der one oscillation is the fraction of the slopes of these 
lines, which turns out to be a constant, irrespectively of 
the method we employ. Since for increasing TV the aver- 
age width of oscillations tends to 0, Ax — > 0, but the area 
under one oscillation remains constant, one might con- 
clude that the height of a wave ought to tend to infinity, 
i.e. Ay — > oo, i.e. the potential has a Dirac-delta-like sin- 
gularity at every point in the TV — > oo limit. In order to 
confirm this conclusion, we have calculated the Fourier- 
transform of sum of sinc-functions centered at random x 
value. The result of this calculation showed qualitatively 
similar behavior as the one presented in Fig. [5] 



In the previous section we demonstrated that a quan- 
tum potential whose energy eigenvalues follow GOE level 
statistics is nowhere differentiable in the TV — > oo limit. 
Can one classify the behavior of this potential any fur- 
ther? Is there a characteristic width of the oscillations 
appearing in the potential? If 'yes', what is the asymp- 
totic value of this characteristic width as TV increases? 



V. CONCLUSIONS 

In this paper we investigated the inverse BGS conjec- 
ture and the validity of the counterexample given by Wu 
et al. The Hamiltonian of our model takes its most stan- 
dard formp 2 /2m+U(x), and we created one-dimensional 
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FIG. 3. (a) The potential calculated with the Wu method for N = 5 (continuous line), N = 50 (dashed line), iV = 250 (points) 
in the interval [0; 2]. (b) The meaning of Ax and Ay. 



• Wu method 




FIG. 4. The average distance of two adjacent extrema, Ax, 
is plotted as a function of N for both methods. The solid 
and dashed lines are plotted only to guide the eye and are 
obtained by fitting the function A + BN~ C onto the data. 
The inset is the same figure on a log-log scale. 



quantum potentials, V(x), whose first N energy levels fol- 
low GOE properties. We investigated the N — > oo limit 
and showed that such a potential is nowhere differentiablc 
and nowhere continuous, moreover may it has a singular- 
ity at "every" point. No sensible physical meaning can be 
associated with such an ill-behaving potential, therefore, 
the question whether it is chaotic or not, cannot be in- 
terpreted either. We mention here, however, the logically 
plausible, but improbable scenario that a more compli- 
cated Hamiltonians, e.g. Harper's model cos (p) +cos (a;) 
[13, HH, may reproduce the GOE statistics. 
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FIG. 5. The average vertical "jump" between two adjacent 
extrema, Ay, is depicted as a function of N obtained using 
both methods. 

We used a method originally employed by Wu et al. , 
and the dressing transformation in order to find a poten- 
tial whose entire spectrum obeys GOE statistics. These 
methods are based on different ideas, nevertheless they 
provided quantitatively the same results for the relevant 
measures investigated. The range of N was between 5 
and 250. 

We further analyzed the N — > oo limit by calculating 
the Fourier spectrum of the fluctuating part of the po- 
tential. We demonstrated that its Fourier components 
do not tend to zero at large magnitudes of /. Therefore 
the potential loses its characteristic length scale, and its 
structure resembles that of a fractal, the dimension is 
above 1.6. This finding also supports our conclusion, 
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namely a potential having GOE spectrum does not have a 
classical analogue, therefore it is not an acceptable coun- 
terexample against the inverse BGS conjecture. 



FIG. 6. (Color online) The histogram of the parameter c is 
depicted for different values of N. 
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FIG. 7. The average of c as a function of TV with the Wu 
method and the dressing transformation. 
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